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We report on the charge kinetics of La 2 Cu0 4 lightly doped with Li and Sr. Impedance spec- 
troscopy measurements down to 25 mK and from 20 Hz to 500 kHz reveal evidence for low energy 
charge dynamics, which slow down with decreasing temperature. Both systems are acutely sensi- 
tive to stoichiometry. In the case of Sr substitution, which at higher carrier concentration evolves 
to a high temperature superconductor, the ground state in the pseudogap-doping regime is one of 
spatially segregated, dynamic charge domains. The charge carriers slow down at substantially lower 
temperatures than their spin counterparts and the dynamics are particularly sensitive to crystallo- 
graphic direction. This is contrasted with the case of Li-doping. 



PACS numbers: 74.72.Dn, 74.25.Fy, 74.25.Nf, 75.50.Lk 
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Hole-doping into the Cu0 2 planes of the antiferro- 
magnetic Mott insulating compound, La 2 Cu0 4 , can be 
achieved for example with Sr or Li. In the case of 
La2Cui_ :r Li E 04 (LLCO), the holes tend to be local- 
ized near the Li-sites [l|, whereas with La2- x Sr a: Cu04 
(LSCO) they are more mobile [2, [3]. In both systems 
doping (a;) suppresses the Neel temperature (Tjv) and a 
short-range "spin-glass"_phase is exp osed near x = 0.03 

at ~ 8K 0, H la 0, 0, S, H, Efl El, ELI El. it is 

only for x > 0.055 that superconductivity emerges in 
LSCO. LLCO, on the other hand, remains an insula- 
tor [12J], yet has a similar magnetic phase diagram at 
low dopings [71, LL3J . The intermediate spin-glass region 
in LSCO allows the study of spin and charge dynam- 
ics prior to the onset of high temperature superconduc- 
tivity (HTS). A wide range of investigations in this in- 
termediate region have provided evidence for the pre- 
sense of antiferromagnetic domains separated by anti- 



phase boundaries [7j, la, LUJ]. In so far as the charge is 



concerned, it remains unclear whether the system is con- 
ducting at these dopings. There have been indications of 
spatial se greg ation of charge into locally ordered regions 

Recent infrared studies are in- 
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consistent with the notion of static charge ordering [21 1. 
These results are supported by low frequency resistance 
noise measurements 22, 23]. 

To elucidate the charge dynamics we performed low 
frequency impedance spectroscopy (20 Hz to 500 kHz) on 
LLCO and LSCO, both with x = 0.03 — a probe sen- 
sitive to correlated domains [24]. Our studies show that 
the charge slows down with decreasing temperature (T) 
in both systems. We observe a dispersion in the con- 
ductivity that we attribute to dielectric losses in both 
systems. We show that LSCO in the pseudogap-doping 
regime studied here manifests spatial segregation of holes 



into locally ordered, low energy, dynamic regions. The 
charge carriers slow at temperatures a couple of orders 
of magnitude lower than their spin counterparts — an 
anomalous situation in glassy systems in general. Signif- 
icantly these low energy charge dynamics occur primarily 
along the Cu0 2 planes; and none of these is evidenced in 
LLCO, which does not superconduct at any doping. 

The LLCO and LSCO (x = 0.03) high-quality single 
crystals were grown using the traveling-solvent floating- 
zone method 13, |25||. In order to eliminate possible hole 
doping by excess oxygen, the crystals were carefully heat 
treated under reducing conditions. The lithium concen- 
trations were estimated to within ±0.003 [13] ■ The crys- 
tal axes were determined by the x-ray Laue backscatter- 
ing technique and samples were cut with a wire saw for 
both in-plane and out-of-plane measurements. Different 
contacts were used to ensure that the results were intrin- 
sic, including evaporated Au and cured Dupont 6838 Ag 
paste. We tried depositing an A1 2 3 layer and subtract- 
ing its contribution, but the results were unchanged. We 
also varied the geometry to confirm that the effects were 
not from depletion layers at the contacts. We performed 
our measurements down to 1.3 K in He dewars and to 
300 mK using an adiabatic demagnetisation refrigerator. 
Measurements were extended to 25 mK using two differ- 
ent dilution refrigerators. We measured the impedance, 
Z and the phase angle, 9. The phase angle varied from 
approximately 0° (resistive) at room temperature and low 
frequency (/ « 100 Hz) to < 90° (low-loss dielectric) at 
low temperatures and higher frequency (/ « 100 kHz). 
The impedance varied by up to nine orders of magni- 
tude in the course of an experiment. To compare such 
extreme variance we extracted from our impedance mea- 
surements the real e' = — ^J' £ n e A and imaginary permit- 



tivity £ 



and the conductivity a = lo£qe'' 



luZe A> 

where d is the distance between the contacts, A is the area 
of the contacts, u) — 2nf, and s Q — 8.854 x 10~ 12 F-m _1 
is the vacuum permittivity [24|. 

Park et al. [26j showed that an electronic glass occurs 
in LLCO (x = 0.023), seen by a steplike drop in e' at 
a characteristic T which increases with /. Having repli- 
cated their results (not shown here), we then studied the 
effect of additional carrier doping. A small increase in 
doping to x — 0.03 causes an astonishing increase in con- 
ductivity and reveals strong peaks in e'(T) (Fig. [1]), the 
temperature of which increases with /. The rounded- 
ness and low frequencies involved suggest polarization of 
electronic domains with a wide range of characteristic 
frequencies. Noise is introduced above the peak T below 
1 kHz, consistent with domain walls switching in random 
steps with the application of a slowly varying field. The 
effect is less noticeable at higher frequencies where the 
domains are excited reversibly if at all or at lower tem- 
peratures where they are sluggish on these time scales. 
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FIG. 1: (Color online) Out-of-plane real permittivity of 
LLCO (x = 0.03) vs T. Inset: T and / dependence of the 
peaks in e' . (1) Arrhenius fit; (2) Vogel-Fulcher and power- 
law fits (see text for details). The error-bars (not shown) are 
smaller than the symbols. 

The / dependence of the peak in e' is fitted against 
three standard forms (Fig. [TJ inset) [24]]. An Arrhe- 
nius fit (exp(— E/UbT)), giving an excitation energy of 
E = 59 K, is more convex than the data points. A Vogel- 
Fulcher fit (exp(—E/kB(T — Tf))) gives a freezing tem- 
perature of Tf = -5.8(7) K < and E = 189(32) K. A 
power-law fit ((T — Tf) n ) gives unphysical parameters 
(n = 7.5). This agrees with the following analysis of the 
conductivity and will be commented on below. 

We extracted the conductivity from our impedance 
measurements (Fig. [21). The linearity of the 20 Hz data 
below *-- 40 K reveals Mott Variable Range Hopping 

d is 




(VRH), a(T) = CT exp(-(To/r)^), where (3 = 



(d+i) 



FIG. 2: (Color online) Out-of-plane o of LLCO vs T" 1/4 . 
In this plot a straight line represents VRH. The extrapolated 
/ independent VRH component we subtract from the data 
is parallel to the double-headed arrow. Lower inset: rep- 
resentative calculated dielectric losses with VRH component 
subtracted (see text for details). Upper inset: T and / de- 
pendence of the peaks in the dielectric losses. The lines are 
the same fits as in Fig. [TJ (inset). The error-bars (not shown) 
are smaller than the symbols. 



the dimensionality of the hopping [27j • There is a / de- 
pendent deviation from VRH. This may be understood in 
terms of two components to the losses: VRH (indicated 
by the double-headed arrow in the main panel of Fig. [2j) 
and dielectric loss. We subtracted the extrapolated best- 
fit / independent VRH component from a to reveal the 
/ dependent loss peaks (Fig. [21 lower inset). A plot of 
the T dependence of these peaks (Fig. [21 upper inset) is 
in reasonable agreement with the data shown in the in- 
set to Fig. [TJ An Arrhenius fit again fails to describe the 
data, and a Vogel-Fulcher fit yields Tf = —6.0(7) K and 
E = 165(22) K. A power law fit again gives unphysical 
parameters (n = 7.7). 

The in-plane measurements (not shown) are qualita- 
tively similar, although a and e' are both an order of 
magnitude larger than the out-of-plane values at low tem- 
peratures. In the limit of / — > our results agree with 
DC measurements on comparable ceramic samples [12J |. 
while at higher frequencies our experiments reveal similar 
low energy charge dynamics to those seen in the out-of- 
plane conductivity above. 

The negative Tf values obtained from the Vogel- 
Fulcher fits suggest that the charge does not freeze in 
these samples, although it continues to slow as T is re- 
duced. It is possible that other fits over a lower temper- 
ature region will yield a different result: this remains to 
be checked with further experiments. Nevertheless, the 



present results differ to those of Park et al. [26fl which 
show Tf of 3 — 5 K. This difference demonstrates the ex- 
treme sensitivity of the coupling between the low energy 
spin and charge dynamics to carrier doping, possibly re- 
lated to the large ^f (x ~ 0.03) [3, |T3| . At x = 0.03 none 
of the standard scaling analyses adequately describe the 
low energy dynamics of the charge over the temperature 
range studied. 
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FIG. 3: (Color online) Real permittivity of LSCO vs /, with 
(a) field in-plane and (b) out-of-plane. The dotted lines are 
guides to the eye. Insets are plots of the T and / dependence 
of the peaks. The line in the lower inset shows two fits (see 
text for details). The error-bars (not shown) are smaller than 
the symbols. 

We now turn to our results on LSCO. Extracting e' 
from the impedance measurements reveals broad T de- 
pendent resonance-like peaks when plotted against / 
(Fig. [3]). The broadening of the dielectric dispersion on 
cooling can be attributed to a broadening of the distri- 
bution of relaxation frequencies. This is suggestive of a 
distribution of sizes of electronic domains acting as co- 
herent resonant oscillators. The peak / decreases when 
T is lowered, indicative of the dynamics slowing down. 
These observations are in accord with the emergence of 
dynamical charge heterogeneities at ultra-low tempera- 
tures. Noise measurements, which reveal distinct switch- 



ing fluctuations at temperatures of < 0.3 K, further con- 
I — I r — I, 

mm 



firm this behavior 

Although the results with field in-plane and out-of- 
plane look qualitatively similar, there are some notable 
differences (Fig. [3]). The in-plane measurements are an 
order of magnitude larger. This may be understood by 
higher charge carrier mobility in the Cu0 2 planes, en- 
hancing the measured polarization of the domains. Sec- 
ondly, the in-plane frequency peaks asymptotically ap- 
proach ~ 1kHz as T — > (Fig. Efa); inset). This sug- 
gests that the sluggishness of the in-plane domains levels 
off as T — > 0. The out-of-plane peaks, on the other hand, 
show a distinct kink at approximately 160 mK (Fig.[3](b); 
inset). We fitted the data below this kink to the three 
standard forms used previously. In this instance a Vogel- 
Fulcher fit reduced to an Arrhenius fit (Tf rs 0), with 
E = 281 mK, and a power-law fit yielded n — 1.56(24) 
and Tf = 42(6) mK. Thus, while the in-plane resonant 
frequency levels off as T — > 0, the out-of-plane dynam- 
ics continue to slow. Whether it has finite freezing at 
Tf = 42 mK (as in the power-law fit) or simply activated 
behavior (as in the Arrhenius fit) remains unclear. 
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FIG. 4: (Color online) a of LSCO vs VRH T exponents, 
(a) P = 1/2, (b) P = 1/3, and (c) /3 = 1/4. Upper curves 
show in-plane data, lower show out-of-plane, which is shifted 
down for clarity. Inset: adjusted in-plane a vs /. The lines 
are guides to the eye. The error-bars (not shown) are smaller 
than the symbols. 

Next we consider the conductivity of this system 
(Fig. 21). When the field is in-plane, there is a simi- 
lar / dependent deviation to that seen in LLCO, but at 
lower temperatures. This lends credence to the low en- 
ergy, dynamic charge heterogeneities mentioned above. 
We note that although the signature for the onset of dy- 
namic charge heterogeneities occurs in the mK regime 
(Fig. [3]), the dispersion shown in Fig. |4] commences at 
temperatures of a few kelvin. We again subtracted the 
VRH component. The losses are found to be relatively 



frequency independent for / < 1 kHz, but highly disper- 
sive for / > 1 kHz (Fig. IU inset). It is possible to dis- 
tinguish a similar but significantly smaller effect in the 
out-of-plane conductivity for T < 500 mK (not visible 
on the scale shown). The deviation occurs first at the 
higher frequencies, and for example at 100 mK is mono- 
tonic with the 320 kHz conductivity 22% larger than the 
120 Hz value. That this effect is relatively small may be 
due to the low out-of-plane charge carrier mobility with 
a concomitant reduction in polarizability. The observed 
anisotropy (Figs [3] and [4]) indicates that the low energy 
dynamics of the electronic domains occur mainly in the 
Cu0 2 planes. 

We then estimated the barrier energy by fitting the 
conductivity to the VRH equation. For the low T in- 
plane region (Fig. (He); T= 1.3-5K&T = 0.04-1K) 
best-fits yield = 1/4 and T = 3.1 x 10 3 K & T = 
340 K. We contrast these x = 0.03 results with x = 0.02 
data from the literature. Ellman et al. [28j reported 
= 1/2 and T Q = 74 K for x = 0.02 at low tempera- 
tures (0.3-10 K), and associated the value with weak 
screening between the interacting charge carriers. To the 
best of our knowledge no low T data oni = 0.02 single 
crystals is available. Nevertheless, this comparison shows 
that a small increase in carrier concentration in this dop- 
ing region is sufficient to cause an onset of screening 
of interactions at low temperatures. This is consistent 
with a rapid doping induced increase in carrier mobil- 
ity [29]. In the lowest temperature out-of-plane region 
(Fig. gta); 0.04-0.12 K) we obtain = 1/2 and T = 7K, 
while at higher temperatures (Fig. [He); 10-50 K) = 1/4 
and Tq = 163 K [30]. Between these temperature ranges 
(Fig.BJb); 0.5-10 K) a best fit yields = 1/3 in agree- 
ment with other measurements in a similar temperature 
range [23j. Again contrasting with the literature, Birge- 
neau et al. J3l| reported = 1/4 and T = 1.1 x 10 6 K 
for x — 0.02 in the temperature region 10-100 K. Thus at 
higher temperatures we see a sharp reduction in barrier 
energy as doping increases from x = 0.02 to x = 0.03. 
We associate the = 1/2 at low T with weak screening 
due to lower carrier mobility in this orientation, and the 
= 1/3 with a gradual transition between the higher 
temperature and low temperature regions. 

In summary, we performed low T, low / impedance 
spectroscopy to compare the charge kinetics of LLCO 
and LSCO (x = 0.03). Both systems exhibit low energy 
charge dynamics, which slow down as T is reduced. Com- 
parison with previous work oni = 0.02 26|, \2Q, [29j, 1 3 lj | . 
shows that both systems undergo a sharp increase in 
mobility as doping increases to X = 0.03. Our study 
of LSCO enabled the identification of a ground state 
with dynamic charge heterogeneities in the doping range 
where the pseudogap and spin- glass phase are present. 
The correlated slow dynamics identified are largely de- 
coupled from their glassy spin counterparts and lie pri- 
marily along the Cu0 2 planes which are responsible for 



the emergence of HTS at higher charge carrier concen- 
trations. Notably this was not so for LLCO. Our ob- 
servations call for a reanalysis of the large region of 
the HTS phase diagram occupied by a spin-glass phase 
(0.01 < x < 0.20) and the role of low energy charge do- 
mains on the emergence of superconductivity in doped 
Mott insulators. 
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